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Abstract 

Di-jet production is studied in collisions of quasi-real photons at e + e~ centre-of-mass ener- 
gies ^/s^ from 189 to 209 GeV at LEP. The data were collected with the OPAL detector. 
Jets are reconstructed using an inclusive fc^-clustering algorithm for all cross-section mea- 
surements presented. A cone jet algorithm is used in addition to study the different struc- 
ture of the jets resulting from either of the algorithms. The inclusive di-jet cross-section is 
measured as a function of the mean transverse energy E 1 ^ of the two leading jets, and as a 
function of the estimated fraction of the photon momentum carried by the parton entering 
the hard sub-process, x y , for different regions of E^ . Angular distributions in di-jet events 
are measured and used to demonstrate the dominance of quark and gluon initiated pro- 
cesses in different regions of phase space. Furthermore the inclusive di-jet cross-section as 
a function of \rf et \ and |Ar/j et | is presented, where rji et is the jet pseudo-rapidity. Different 
regions of the x^-x~-space are explored to study and control the influence of an underlying 
event. The results are compared to next-to-leading order perturbative QCD calculations 
and to the predictions of the leading order Monte Carlo generator PYTHIA. 
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1 Introduction 



We have studied the production of di-jets in the collisions of two quasi-real photons at an 
e + e~ centre-of-mass energy ^fs^, from 189 to 209 GeV, with a total integrated luminosity 
of 593 ptT 1 collected by the OPAL detector at LEP. Di-jet events are of particular interest, 
as the two jets can be used to estimate the fraction of the photon momentum participating 
in the hard interaction, which is a sensitive probe of the structure of the photon. The 
transverse energy of the jets provides a hard scale that allows such processes to be calculated 
in perturbative QCD. Fixed order calculations at next-to-leading order (NLO) in the strong 
coupling constant a s for di-jet production are available and are compared to the data, 
providing tests of the theory. Leading order Monte Carlo (MC) generators are used to 
estimate the importance of soft processes not included in the NLO calculation. 

Inclusive jet cross-sections in photon-photon collisions have previously been measured 
at v/s^ = 58 GeV at TRISTAN Pig and at ^ from 130 to 172 GeV at LEP |3E|. This 
paper extends the latter analysis to higher e + e~ centre-of-mass energies, and provides an 
approximately thirty-fold increase in integrated luminosity. The /c_i_-clustering algorithm [S] 
is used as opposed to the cone algorithm [U] in [3J E] for the measurement of the differential 
cross-sections, because of the advantages of this algorithm in comparing to theoretical 
calculations [7j. The cone jet algorithm is used to demonstrate the different structure of the 
cone jets compared to jets defined by the /cj_-clustering algorithm. The large amount of data 
allows us to measure the cross-section for di-jet production in photon-photon interactions 
as a function of the mean transverse jet energy E^, the jet pseudo-rapidity |^ et | and the 
absolute difference in pseudo-rapidity lA^ 6 *! of the jets, with rf et = — lntan(^ et /2) . For 
the first time, the differential cross-section is also measured as a function of the estimated 
fraction of the photon momentum carried by the parton entering the hard sub-process, x 7 , 
with full unfolding for detector effects. Angular distributions in di-jet events are measured 
and used to demonstrate the dominance of quark and gluon initiated processes in different 
regions of phase space. 

At e + e~ colliders the photons are emitted by the beam electrons 2 . Most of these photons 
carry only a small negative four-momentum squared, Q 2 , and can be considered quasi-real 
(Q 2 ~ 0). The electrons are hence scattered with very small polar angles and are not 
detected. Events where one or both scattered electrons are detected are not considered 
in the present analysis. Three processes contribute to di-jet production in photon-photon 
collisions: the direct process where two bare photons interact, the single- resolved process 
where a bare photon picks out a parton (quark or gluon) of the other photon, and the 
double-resolved process where partons of both photons interact This separation is 
only unambiguous at leading order. At higher orders it becomes dependent on the process 
scales. 

1 The coordinate system of OPAL has the z axis along the electron beam direction, the y axis pointing 
upwards and x towards the center of the LEP ring. The polar angle 9 and the azimuthal angle cf> are 
defined relative to the +z-axis and +x-axis, respectively. 

2 Positrons are also referred to as electrons. 
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2 The OPAL detector 



A detailed description of the OPAL detector can be found in Only the main features 
relevant to the present analysis will be given here. 

The central tracking system is located inside a solenoidal magnet which provides a 
uniform axial magnetic field of 0.435 T along the beam axis. The magnet is surrounded in 
the barrel region (| cos#| < 0.82) by a lead glass electromagnetic calorimeter (ECAL) and a 
hadronic sampling calorimeter (HCAL). Outside the HCAL, the detector is surrounded by 
muon chambers. There are similar layers of detectors in the endcaps (0.82 < | cos 6\ < 0.98). 
The small angle region from 47 to 140 mrad around the beam pipe on both sides of the 
interaction point is covered by the forward calorimeters (FD) and the region from 33 to 
59 mrad by the silicon tungsten luminometers (SW). 

Starting with the innermost components, the tracking system consists of a high precision 
silicon micro-vertex detector, a vertex drift chamber, a large volume jet chamber with 159 
layers of axial anode wires and a set of z chambers measuring the track coordinates along 
the beam direction. 

The barrel and endcap sections of the ECAL are both constructed from lead glass 
blocks with a depth of 24.6 radiation lengths in the barrel region and more than 22 radiation 
lengths in the endcaps. The HCAL consists of streamer tubes and thin multiwire chambers 
instrumenting the gaps in the iron yoke of the magnet, which provides the absorber material 
of 4 or more interaction lengths. 

The FD consists of two cylindrical lead-scintillator calorimeters with a depth of 24 
radiation lengths divided azimuthally into 16 segments. The SW detectors consist of 19 
layers of silicon detectors and 18 layers of tungsten, corresponding to a total of 22 radiation 
lengths. 



3 Monte Carlo simulation 

The MC generators PYTHIA 5.722 PHE] and PHOJET 1.10 [H] are used to study de- 
tector effects. PYTHIA is based on leading order (LO) QCD matrix elements for massless 
quarks with added parton showers and hadronisation. PHOJET also simulates hard in- 
teractions through perturbative QCD in LO, but includes soft interactions through Regge 
phenomenology before the partons are hadronised. The probability of finding a parton 
in the photon is taken from parametrisations of the parton distribution functions. The 
default choices of SaS ID [T3| for PYTHIA and LO GRV [H] for PHOJET are taken for 
the samples used to study detector effects. 

An increased flow of transverse energy, E^, apparently not directly related to the hard 
subprocess has been observed in photon-hadron scattering [T^j, and has been labeled the 
underlying event. Both PHOJET and PYTHIA include a model of multiple parton inter- 
actions (MIA) to simulate such effects. In PYTHIA the amount of MIA added to the event 
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is controlled by a lower cutoff parameter p™ 1 , which describes the transverse momentum of 
the parton involved. Following the studies carried out in [31 , P™ 1 is set to 1.4 GeV for 
the SaS ID parton densities. In PHO JET the default setting for MIA is used. 

Three non-signal processes are important: hadronic decays of the Z°, where initial state 
photon radiation has reduced the centre-of-mass energy of the hadronic final state to be 
close to the Z° mass, 77 — > tt reactions, and photon-photon collisions where one of the 
photons is virtual (7*7) but the scattered electron is not detected. The hadronic Z° decays 
are simulated using PYTHIA. The pair-production of r-leptons in photon-photon collisions 
is simulated using BDK [16J. Deep-inelastic electron-photon scattering events are studied 
with the HERWIG 5.9 [Tlj generator. 

All signal and background MC samples were generated with the full simulation of the 
OPAL detector ^H]- They are analysed using the same reconstruction algorithms as are 
applied to the data. 



4 Definition of di-jet observables 

All cross-section measurements use jets reconstructed with the inclusive fcj_-clustering al- 
gorithm as proposed in with Rq = 1. In addition, a cone jet algorithm jB] with a cone 
size of 1.0 in ?7-0-space is employed to study the dependence of the jet structure on the 
algorithm used. A di-jet event is defined as an event with at least two jets fulfilling the 
requirements detailed below. In events with more than two jets, only the two jets with the 
highest values are taken. 

The primary intentions of this analysis are to study the ability of QCD theory to 
describe jet production in photon-photon collisions, and to explore the photon structure 
revealed in these hadronic interactions. The most advanced theoretical predictions to 
date are provided by fixed order perturbative calculations up to NLO for the production 
of di-jets. These calculations need as input a scale that is in principle arbitrary, but 
commonly set to a value related to the hardness of the interaction. Possible choices for 
di-jet production are for example the mean transverse energy or the maximum E^ of 
the di-jet system. 

The separation of quasi-real and virtual photons is somewhat arbitrary and therefore 
needs to be defined. For this analysis we choose values of Q 2 < 4.5 GeV 2 to define quasi- 
real photons. It is the same value as used in previous analyses [3J El , motivated by the 
acceptance of the low angle calorimeters. The median Q 2 resulting from this definition 
cannot be determined with the data since the scattered electrons are not tagged. For the 
kinematic range of this analysis both PHOJET and PYTHIA predict the median Q 2 to be 
of the order 10" 4 GeV 2 . 
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4.1 Properties of di-jet events 



In LO QCD, neglecting multiple parton interactions, two hard parton jets are produced 
in 77 interactions. In single- or double-resolved interactions, these jets are expected to be 
accompanied by one or two remnant jets. A pair of variables, x+ and x~ can be defined [TH] 
that estimate the fraction of the photon's momentum participating in the hard scattering: 

E (£ jct +^ ct ) E (£ jct -^ ct ) 

I jets=l,2 . _ jets=l,2 , , 

xZ = — and x. ~ 



7 ~ J2(e+ Pz ) 7 ^£(E- Pz ) ' 

hfs hfs 

where p z is the momentum component along the z axis of the detector and E is the 
energy of the jets or objects of the hadronic final state (hfs). In LO, for direct events, 
all energy of the event is contained in two jets, i.e., x+ = 1 and x~ = 1, whereas for 
single-resolved or double-resolved events one or both values are smaller than 1. The di-jet 
differential cross-section as a function of x 7 is therefore particularly well suited to study 
the structure of the photon, since it separates predominantly direct events at high x 7 
(x 7 > 0.75) from predominantly resolved events at low a; 7 (x 7 < 0.75). The fraction of 
direct, single-resolved and double-resolved events as a function of x 7 predicted by PYTHIA 
is shown in Figure Q (a)-(c). The dominance of resolved events for x 7 < 0.75 is clearly 
visible. In these distributions and in the definitions below, x 7 indicates that each event 
enters the distribution twice, at the value of 2+ and the value of x~. 

Due to the different nature of the underlying partonic process one expects different 
distributions of the angle G* between the jet axis and the axis of the incoming partons 
or direct photons in the di-jet centre-of-mass frame. The leading order direct process 
77 — > qq proceeds via the t-channel exchange of a spin-| quark, which leads to an angular 
dependence oc (1 — cos 2 ©*) . In double resolved processes the sum of all matrix elements, 
including a large contribution from spin-1 gluon exchange, leads to an approximate angu- 
lar dependence oc (1 — |cos©*|) -2 [2U|. The contribution of the different processes to all 
resolved events depends on the parton distribution functions of the photon. An estimator 
of the angle 9* can be formed from the pseudo-rapidities of the two jets as 

cos6* = tanh \ i j , (2) 

where it is assumed that the jets are collinear in <ft and have equal transverse energy. 
Only |cos©*| can be measured, as the ordering of the jets in the detector is arbitrary. 
To obtain an unbiased distribution of I cos 0*1 the measurement needs to be restricted to 



the region where the di-jet invariant mass My = 2E^ t /yl — |cosO*| 2 is not influenced 

by the cuts on E 3 ^ 1 In the present analysis a cut of Mjj > 15 GeV ensures that the 
|cosO*| distribution is not biased by the restrictions on E^ for the range |cos©*|<0.8 and 

\ff et \ = | ^j 6 * + ^2 6t ) /2| < 1 confines the measurement to the region where the detector 

resolution on |cosO*| is good. 
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4.2 Differential di-jet cross-sections 



The following differential cross-sections are measured, where the labels 1 and 2 refer to the 
two jets with highest Ei£ in the event, defined by the k± algorithm: 

^ with 4 et = ^ + E ^ 2 and Elf > 5 GeV (3) 
dE^ T 2 1 ; 



dfTrfii 



dijct 



dx 7 



dfTHii 



dijct 



dlogio (x 7 ) 

dcT di j et dddijet dddijet 



diu&r dic^r d|A^i 

dcr di j et 



dlcos 0*1 



in 3 bins of E^ [5 - 7 - f f - 25] GeV (4) 

for 5 GeV < E 3 * < 7 GeV (5) 

for E$* > 5 GeV (6) 

for E 3 * > 5 GeV, |r/ jct | < 1, M$ > 15 GeV (7) 
with in all cases 



et \ E t\ ~ E t\\ 1 

|r7i C 9 | < 2 and — Hr ^— < - 

1/1,21 4 



Here, \rf^ ntr \ and |^ d | denote the jet with the smaller and larger value of \r]^ et \ respectively, 
and | Ar/j et | is defined to be the absolute distance in pseudo-rapidity between the two leading 
jets. 

The combination of the second condition in Equation (jHJ) with the minimum E 3 ^ re- 
quirement defines asymmetric E 3 ^ thresholds for the two jets of the di-jet system, which 
is important in comparisons to NLO QCD calculations [21 j . This method of defining 
asymmetric thresholds has previously been used in [2*2] . 

Four regions in x^-x'-space are considered (see Figured] (d)): (A) the complete 
space (full x^ range), (B) both x+ and x~ larger than 0.75 (x^ > 0.75), (C) either x+ or 
x~ smaller than 0.75 (x* or x~ < 0.75), (D) both x^ and x~ smaller than 0.75 (xt < 0.75). 

The cross-sections (J3j), (jlj) and (jSJ) are measured in regions (A), (C) and (D). For the 
cross-sections in (jBJ regions (C) and (D) are considered. The cross-section as a function of 
|cos0*| in (J2J) is measured in regions (B) and (D). 



4.3 Jet structure in di-jet events 

The internal structure of jets is studied using the jet shape, which is defined as the fractional 
transverse jet energy contained in a subcone of radius r concentric with the jet axis, 
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averaged over all jets of the event sample: 



1 7^-1 jet / \ 

^ r )^^E F jeV ^ With r=^{A^T{W 2 - (9) 
N i«* E t ( r = 1-0) 

iVjet is the total number of jets analysed. Both k± and cone jets are analysed in this way. 
As proposed in [221, only particles assigned to the jet by the jet finders are considered. 
Events entering the jet shape distributions are required to have at least two jets with a 
transverse energy 3 GeV < E^ < 20 GeV and a pseudo-rapidity \r]^ ct \ < 2. The cone jet 
algorithm is used in addition to the /cj_-clustering algorithm to demonstrate the different 
structure of the cone jets with respect to those defined by the fc^-clustering algorithm. 

The jet shape is measured in the two regions of x+-x~- space, x^ < 0.75 and x^ > 0.75, 
in four bins of with bin boundaries at 3, 6, 9, 12 and 20 GeV and four bins of \r]^ et \ 
between and 2, each bin 0.5 units wide. 



5 Event selection 



In this analysis, a sum over all particles in the event or in a jet means a sum over two 
kinds of objects: tracks and calorimeter clusters, including the FD and SW calorimeters. 
A track is required to have a minimum transverse momentum of 120 MeV and at least 
20 hits in the central jet chamber. The point of closest approach to the origin must have 
a distance of less than 25 cm in z and a radial distance of less than 2 cm to the z-axis. 
Calorimeter clusters have to pass an energy threshold of 100 MeV in the barrel section or 
250 MeV in the endcap section for the ECAL, 600 MeV for the barrel and endcap sections 
of the HCAL, 1 GeV for the FD, and 2 GeV for the SW. An algorithm is applied to 
avoid double- counting of particle momenta in the central tracking system and their energy 
deposits in the calorimeters j^j. The measured hadronic final state for each event consists 
of all objects thus defined. 

Di-jet events are preselected using the k±_ algorithm by requiring at least two jets 
with |^ et | < 2 and a transverse energy E^ > 3 GeV. Photon-photon scattering events are 
selected using the requirements detailed below. The corresponding distributions in Figured 
compare the sum of the simulated signal and background processes to the data, uncorrected 
for detector effects. For each distribution shown, all selection criteria are applied except the 
one on the quantity plotted. The signal MC generators PHOJET and PYTHIA are found 
to underestimate the cross-section by about 20% in these comparisons, and are scaled up 
accordingly. Of all non-signal processes studied, only those listed in Section |3] contribute 
significantly. Comparisons of the rate of di-jet events in photon-photon collisions where 
one of the photons is virtual (see for example [21]) show that the prediction of the MC 
generator used is too low by about a factor of two. The prediction of the contribution from 
7*7 events has been scaled up accordingly. 

All distributions are sufficiently well described by the sum of signal and background 
contributions. The total contribution of non-signal processes to the selected event sample 
is about 5% after the following selection criteria have been applied: 
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• The sum of all energy deposits in the ECAL and HCAL (Figure 121(a)) has to be 
less than 55 GeV to remove background from hadronic Z°-decays in events with a 
radiative return to the Z°-peak. 

• The visible invariant mass measured in the ECAL, Wecal, has to be greater than 
3 GeV to suppress low energy events. 

• The missing transverse momentum of the event, Pt,miss ; calculated from the mea- 
sured hadronic final state, has to be less than 0.05 • -Ebeam- 

• At least 7 tracks must have been found in the tracking chambers. This cut re- 
duces mostly the contamination from 77 — > tt events. The distribution of the track 
multiplicity is shown in Figure 121(b). The discrepancy in shape between data and 
simulation is not present when using PYTHIA instead of PHOJET as signal MC 
generator, and is addressed in the study of model dependencies in Section |H1 

• To remove events with scattered electrons in the FD or in the SW calorimeters, 
the total sum of the energy measured in the FD has to be less than 0.25 ■ -Eream 
and the total sum of the energy measured in the SW calorimeter has to be less than 
0.18 • .Ebeam- These cuts also reduce the contamination from hadronic Z°-decays with 
their thrust axis close to the beam direction. The energy sum in the FD calorimeter 
scaled by the beam energy is shown in Figure 121(c). 

• The z position of the primary vertex is required to satisfy \z\ < 5 cm and the net 
charge Q of the event calculated from adding the charges of all tracks is required to 
be \Q\ < 5 to reduce background due to beam-gas interactions. 

• To remove events originating from interactions between beam electrons and the beam- 
pipe the radial distance of the primary vertex from the beam axis has to be less than 
3 cm. 

• To further reject background from hadronic Z°-decays and from deep-inelastic electron- 
photon scattering an invariant mass, Mjih2, is calculated from the jet with highest 

in the event and the four-vector sum of all hadronic final state objects in the 
hemisphere opposite to the direction defined by this jet. The quantity Mjih2 is a 
simple reconstruction of the Z°-mass in case of background from hadronic Z°-decays, 
and will therefore be larger on average for this type of background than for signal 
events. Events with Mjih2 > 55 GeV are rejected. The distribution of Mjim is 
shown in Figure 121(d). 

We use data at centre-of-mass energies y/s^ from 189 GeV to 209 GeV. For the purpose 
of this analysis, the difference between the data taken at the various values of ^/s^ is small 
and therefore the distributions for all energies have been added. The luminosity weighted 
average centre-of-mass energy y/s^ is approximately 198.5 GeV. The efficiency to trigger 
di-jet events in the region of phase space considered in this analysis has been shown to be 
close to 100% H- 
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6 Transverse energy flow in di-jet events 



NLO QCD calculations do not take into account the possibility of an underlying event 
which leads to an increased .Bp-flow and therefore to an increased jet cross-section above 
a given threshold of E^ . In PYTHIA and PHOJET the underlying event is simulated 
by multiple parton interactions. The contribution from multiple parton interactions is not 
known a priori, but has to be adjusted to give a good description of the data. In this 
analysis the size of this contribution is taken from our previous study of di-jet events in jlj. 
The transverse energy flow from an underlying event is expected to be small compared to 
the transverse energy of the leading jets, and it is not correlated to the direction of the 
jet axes. The energy flow outside the jets will therefore be most sensitive to the presence 
of an underlying event |22 3 - Additional energy outside the leading jets will shift the x 7 
distributions towards lower values. 

To study the performance of the MC generators in describing the energy flow several 
uncorrected distributions are used. The average B T -flow per event is measured with respect 
to the jet function of A0 and A77. The variable 17 is equivalent to 77, except that it 

is signed positively if 2+ is greater than or equal x~ and signed negatively otherwise. The 
definition of fj ensures that the energy flow associated with the "more resolved" photon, 
i.e., the smaller value of x 7 , will always appear on the left hand side of the plots. The 
profiles in A0 consider a range of \Arj\ = 1 around the jet-axis, while a |A0|-range of 7r/2 
around the jet-axis is considered for the profiles in Af). The two leading jets in E^ 1 in each 
event are considered. Another sensitive variable is the energy flow with the leading two 
jets in the event removed, E out , as a function of fj . All objects are excluded inside a cone 
of radius 1.3 in rj-cj) around the two leading jets. 

Multiple parton interactions are important in interactions where the photon is resolved. 
It is therefore interesting to study separately the three cases of (a) two resolved photons, 
(b) one resolved photon, or (c) no resolved photon in the interaction. Experimentally these 
situations can be approximated by choosing events with x^ < 0.75, x^ or x~ < 0.75, or 
x^ > 0.75, i.e. regions (D), (C) and (B) defined in Section HJ 

Figures El and 0] show the jet profiles and the energy flow outside the leading jets. The 
data are compared to a mixture of signal (PHOJET or PYTHIA) and background MC 
simulation. The contributions of signal and background are weighted according to their 
cross-section in each region of phase space. The background MC generators used are the 
same as in Figure 121 

For the 0-profiles in Figure El it is evident that both PHOJET and PYTHIA are able 
to describe the data in the region of the high Et jets around zero. Moving away from 
the jet axis PHOJET predicts an energy flow which is too low compared to the data, 
especially for x^ < 0.75. This corresponds to the area where effects from an underlying 
event are expected to be most prominent. PHOJET improves towards higher x T PYTHIA 
reproduces the data reasonably well. Such differences between PHOJET and PYTHIA are 
not evident in the ^-profiles. Here in the case of resolved photons the E^-How is dominated 
by the photon remnant (s), and is reasonably well described by both generators. The jets 
entering Figure El are selected from the range 10 < E^ < 25 GeV. No significant deviation 



11 



from the behaviour just described is observed when selecting < 10 GeV. 

The energy-flow outside the two leading jets is shown in Figure 0] Again PYTHIA 
describes the data well, while PHOJET is somewhat low. With both models used to 
unfold the data as a systematic check, we conclude that the details of the energy flow 
around the two leading jets are sufficiently well under control and remaining influences are 
included in the systematic uncertainty of the cross-section and jet shape measurements. 



7 Data corrections 

An example of the uncorrected distributions as a function of x 7 including the contribution 
of the remaining background events is shown in Figure 03 To obtain jet cross-sections 
which can be compared to theoretical calculations, we use MC simulations to correct for 
the selection cuts, the resolution effects of the detector and the background from non- 
signal processes. Backgrounds are first subtracted bin-by-bin from all distributions. For 
the differential cross-sections as a function of and x 1 , sizable migration and resolution 
effects are to be expected. We therefore apply a matrix unfolding method, as implemented 
in the GURU program 26J, for these distributions. The results are cross checked using 
a bin-by-bin correction. By definition X<y Cclll only take values between zero and unity. 
At either extremity no adjacent bins are available. To avoid instabilities expected from 
the smoothing procedures in the unfolding for the first and last bin, the central values 
for these bins are taken from the bin-by-bin correction method. The |cosG*|, \rj^ et \ and 
|Ar/j ct | distributions are corrected bin- by-bin, as only small migrations are expected here. 
PHOJET is used as the default signal MC generator for the unfolding. 

The correction method employed for the jet shapes is a bin-by-bin correction using the 
MC simulations to correct for detector effects. The contribution of the same background 
processes as for the cross-section measurements was studied. The influence on the signal 
was found to be less than 1%. Therefore the subtraction of the background was omitted 
in this analysis. Both PYTHIA and PHOJET were used to estimate the correction factors 
to study their model dependence. 



8 Systematic uncertainties 

The overall systematic uncertainty is determined from the sources listed below added in 
quadrature. The same sources are considered for the measurement of the differential cross- 
sections and the jet shapes, with the exception of the background, which has been neglected 
for the jet shapes as discussed in Section 

• To assess the uncertainty associated with the subtraction of background events, the 
predictions for hadronic decays of the Z° and for 77 — > tt reactions are conserva- 
tively varied by 10% without contributing significantly to the systematic error. The 
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prediction of the contribution from 7*7 events has been scaled up by a factor of two 
as described above. By comparing the predictions to the data for large Efd/E-qeam 
and Mj!H2 (see Figure |2J), where this background dominates, we determine that this 
scaling factor can be varied by no more than about 30% in order to keep a good 
description of the data. The scaling factor is varied accordingly. The uncertainty 
from all the background subtraction is typically 2-4%. 

• To estimate the systematic error arising from the specific model used for the un- 
folding, both PYTHIA and PHOJET are used to unfold the data. The estimated 
uncertainty derived from this study is typically 10%, and up to 20% in some cases 
for the differential cross-sections, and 1-2% for the jet shapes. 

• The absolute energy scale of the ECAL calorimeter is known to about 3% [21] for the 
range of jet energies in this analysis. To estimate the influence on the observables the 
energy scale is varied by this amount and the analysis is repeated. The cross-sections 
change by 5-10% due to this variation. The estimated uncertainty for the jet shapes 
is about 1%. 

• The selection criteria described in SectionEJare varied simultaneously both to be more 
restrictive and to allow more events into the analysis to exclude a strong dependence 
on the event selection. Selection criteria based on energy measurements are varied by 
10% of their central value, which is considered conservative given the uncertainty in 
the energy scale and the energy resolution of the calorimeter. The number of tracks 
required and the maximum net charge of the event are changed by ±1. The allowed 
radial distance and z position of the primary vertex are varied by 0.5 cm and 1 cm 
respectively. The uncertainty on the cross-section derived from all these variations is 
typically 5-10%, and up to 20% in some cases for the differential cross-section, and 
about 2-4% for the jet shapes. 

The uncertainty on the determination of the integrated luminosity is much less than 1%, 
and is neglected. For the differential cross-sections the systematic uncertainties evaluated 
for each bin were averaged with the results from its two neighbours (single neighbour for 
endpoints) to reduce the effect of bin-to-bin fluctuations. 



9 Hadronisation corrections 

The differential di-jet cross-sections measured are compared to NLO QCD calculations 
which predict jet cross-sections for partons, whereas the experimental jet cross-sections are 
presented for hadrons. Effects due to the modelling of the hadronisation process are not 
taken into account in the NLO calculation. Because the partons in the MC generators and 
the partons in the NLO calculations are defined in different ways there is as yet no rigorous 
procedure to use the MC generators to correct the data so that they can be compared to the 
NLO parton level predictions. However, as the MC generators are the only available option 
so far, they are used to study the approximate size of these hadronisation corrections. For 
this purpose the prediction of the MC generators at the level of the partonic final state 



13 



is calculated and divided by the prediction obtained from the hadronic final state. The 
resulting correction factor is labeled (1 + Shadr)- The partonic final state consists of all 
partons at the end of the parton shower. The hadronic final state utilises all charged and 
neutral particles with lifetimes greater than 3 x 10 -10 s, which are treated as stable. 

Examples of hadronisation corrections estimated by PYTHIA 6.161 and HERWIG 6.1 
for the observables defined in Section H] are shown in Figure E3 The numerical values can 
be found in j2H|- In PYTHIA the partonic final state is hadronised according to the string 
fragmentation model, while HERWIG uses cluster fragmentation. In all plots the full x+- 
x~-range is considered. The theoretical calculations are corrected bin-by-bin using the 
mean of the hadronisation corrections estimated using PYTHIA and HERWIG. 

Figure |H1 (a) shows (1 + Shadr) as a function of . The correction is less than 10% for 
greater than about 10 GeV, but rises to about 25% for PYTHIA and 15% for HERWIG 
towards small E^ . The corrections for observbles involving the jet pseudo- rapidities are 
dominated by the low E^ region. They are essentially flat and around 20% for PYTHIA. 
HERWIG estimates these corrections to be around 10%. 

Figure El (b) shows (1 + Shadr) as a function of x 1 for the lowest bin in E^ defined 
in Section From the figure it is evident that hadronisation causes large corrections for 
a; 7 > 0.75. The effect is reduced for higher values of E^ , as can be seen in Figured (c) and 
(d). The large corrections for x 7 > 0.75 are mainly due to the large influence hadronisation 
has on the distribution of direct events, which are peaked at x 7 = 1 for the partonic 
final state of the LO calculation, but are much more smeared out at the level of stable 
hadrons (see Figure ^ (a)-(c)). While both the measurement and the NLO calculation are 
perfectly valid for the presented bin sizes, the hadronisation corrections needed for the 
comparison introduces a large migration between the two bins above x 7 = 0.75. For a 
sensible comparison in this region one should therefore consider the sum of the two bins 
above x 7 = 0.75 and compare it to the corresponding sum for the data. 



10 Results 

10.1 Jet structure in di-jet events 

In Figure |7|(a) the jet shape, ^(r), is shown for the k± algorithm for both x^ > 0.75 and 
x^ < 0.75. The first sample is dominated by direct photon-photon interactions and hence 
by quark-initiated jets. As is demonstrated in the figure, jets in this sample are more 
collimated than for small values of x^, where the cross-section is dominated by resolved 
processes and hence has a large contribution from gluon-initiated jets. In both cases the 
jets become more collimated with increasing transverse energy, as is shown in Figure [71(c) . 
There is no significant dependence on the jet pseudo-rapidity (Figure[7|(d)). Both PHOJET 
and PYTHIA give an adequate description of the jet shapes as can be seen in Figures |7|(b), 
(c), and (d). 

Figure[S]compares the shapes of jets defined by the cone algorithm and the k± algorithm, 
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in each case compared to the shape as obtained from PYTHIA. As for the /c_i_-jets, the 
jets defined by the cone algorithm are more collimated in the quark-dominated sample 
and always become more collimated for increasing transverse energy, while there is no 
dependence on the jet pseudo-rapidity The cone-jets are significantly broader than the jets 
defined by the k± algorithm at low E^ . With increasing E^ , jets become more collimated 
and the two jet algorithms give similar results. While the fcj_-jets are well described by 
PYTHIA and PHOJET, the jet shapes obtained for the cone-jets are somewhat broader 
than in the data. 



10.2 Differential di-jet cross-sections 

Only the k± jet algorithm is used for the measurement of the differential di-jet cross- 
sections. The experimental results are compared to a perturbative QCD calculation at 
NLO [29J which uses the GRVHO parametrisation of the parton distribution functions of 
the photon , and was repeated for the kinematic conditions of the present analysis. The 
renormalisation and factorisation scales are set to the maximum E^ in the event. The 
calculation was performed in the MS-scheme with five light flavours and Aqq D = 130 MeV. 
The average of the hadronisation corrections estimated by PYTHIA and HERWIG have 
been applied to the calculation for this comparison. 

Figure Inland Tabled show the differential di-jet cross-section as a function of |cosG*| 
for both x^ > 0.75 and x^ < 0.75. The steeper rise with increasing |cosG*| from the 
dominating spin-1 gluon exchange in the second sample is clearly visible (see Section 0]). 
The shape of both samples is well described by NLO QCD. For x^ < 0.75 the NLO 
calculation is about 20 % below the data. It should be noted that in this region the 
contribution from the underlying event, not included in the calculation, is expected to be 
largest, as discussed in more detail below. For x^ > 0.75 the NLO QCD prediction is about 
20 % above the data. While here the contribution from MIA is small, this region is affected 
by rather large hadronisation corrections as discussed in Section U which translates into 
an uncertainty of the normalisation in comparing the theoretical prediction to the data. 

The differential di-jet cross-section as a function of the mean transverse energy E^ 
of the di-jet system is shown in Figure EH and Table H At high E^ the cross-section 
is expected to be dominated by direct processes, associated with the region x^ > 0.75. 
Consequently we observe a significantly softer spectrum for the case x 7 < 0.75 than for 
the full The calculation is in good agreement with the data for the full x 7 - 

x~-range and for or x~ < 0.75. The cross-section predicted for x 7 < 0.75 is again below 
the measurement. PYTHIA 6.161 is in good agreement with the measured distributions 
using the SaS ID parton densities. 

Figure ^2 and Tables E2 and |H1 show the di-jet cross- section as a function of x 1 and 
logi (x 7 ). The cross-section for the lowest values of E^ shows the largest fraction of 
events at x 7 < 0.75 of the three ranges considered, and is hence most sensitive to gluon- 
initiated processes. The di-jet cross-section logarithmic in x 7 emphasises the region of 
lowest accessible x 7 , which extends down to approximately 0.02. 
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As E 3 ^ increases, the fraction of events with x 1 > 0.75 increases. In MC simulations 
these are predominantly direct events. The sensitivity to the gluon density in the photon is 
hence expected to decrease with increasing E 3 ^ . On the other hand, NLO QCD predictions 
which use E 3 ^ as the process-relevant scale are expected to become more reliable as this 
scale increases. It is hence important to provide measurements at both low and high values 
of -Rp > to study all aspects of the theory. 

PYTHIA using SaS ID is in good agreement with the measured distributions, with a 
tendency to be too low for small values of x 7 . The shaded histogram at the bottom of each 
plot indicates the MIA contribution to the PYTHIA prediction. The numerical values of 
this contibution can be found in [2H]- NLO QCD predicts the shape of the cross-sections 
well for x 7 < 0.75, but is too low by about 10-20% especially at low Eif?. As MIA are 
not included in this calculation it is interesting to note that the MIA contribution to the 
cross-section as obtained from PYTHIA is similar in size to the discrepancy. 

The region of x 1 > 0.75 suffers from large hadronisation corrections as discussed in 
Section El The uncertainty for the data-theory comparison associated with these large 
corrections can be reduced by considering the sum of the two bins above x 1 = 0.75, for 
which NLO QCD indeed gives an adequate description of the data. 

Figure ED demonstrates the effect of using different parton distribution functions of the 
photon on the NLO QCD prediction. AFG HO |S| and GS96 HO |3TJ are used in addition 
to the default GRVHO. The sensitivity of the cross-section to the different gluon density 
in each case is clearly visible for the gg-contribution (Figure H21 (b)), but is less pronounced 
for the full cross-section as can be seen in Figure IT2l (&), due to compensating effects from 
processes involving the quark densities (Figures H21 (c) and (d)). A global analysis, beyond 
the scope of this paper, using for example measurements to constrain simultaneously 
the quark densities hence promises to yield the highest sensitivity to the gluon density in 
the photon. 

In Figure IT31 (Tables El and the same cross-sections as in Figure 1TT1 are shown for the 
case 2+ or x~ < 0.75. Here the cross-section is dominated by interactions where one of 
the two incident photons is resolved. The multiple parton interactions used in PYTHIA 
to model an underlying event are much suppressed in this case, and NLO QCD describes 
both the shape and normalisation of the data well. The opposite effect can be observed 
in Figure ITU (Tables 151 and . where for the case x^ < 0.75 one expects a large influence 
of multiple parton interactions, as demonstrated again by the shaded histogram at the 
bottom of each plot. The cross-sections change by as much as 50% in the low E 3 ^ region, 
when MIA are switched on. For higher E 3 ^ the influence is not as strong. Even with MIA 
switched on, PYTHIA using SaS ID is too low. The deficit visible in the normalisation of 
the NLO calculation is again of similar size as the MIA contribution to the cross-section 
obtained from PYTHIA. 

Complementary information can be obtained by measuring the angular distributions of 
the two highest E 3 ^ jets in di-jet events. The I?? 1 ' 6 *! and lA^ 6 *! dependence of the di-jet 
cross-section is dominated by the low E 3 ^ events. The cross-sections measured are listed 
in Tables |H1 and El In Figure ITS! the di-jet cross-sections as a function of |Ar/j ct |, \rf^ tT \ 
and \ff^ A \ are shown for the case x^ or x~ < 0.75. Again the multiple parton interactions 
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used in PYTHIA to model an underlying event are much suppressed in this case. PYTHIA 
using SaS ID is about 20% too low. The prediction of NLO QCD is in good agreement 
with the data in both shape and normalisation. 

In Figure EH the same cross-sections are presented for the case of x^ < 0.75. As 
expected, the effect of including MIA in PYTHIA is again sizable. When multiple parton 
interactions are switched on, the prediction obtained from PYTHIA reproduces the data 
reasonably well. Again the prediction of NLO QCD is too low by about the size of the 
MIA contribution to the cross-section obtained by PYTHIA. 



11 Conclusions 

We have studied di-jet production in photon-photon interactions with the OPAL detector 
at e + e~ centre-of-mass energies t/s^ from 189 to 209 GeV with an integrated luminosity 
of 593 pb _1 . The data are combined into one sample with a luminosity weighted average 
centre-of-mass energy of approximately y/s^, = 198.5 GeV. Jets are reconstructed using an 
inclusive /^-clustering algorithm for the measurement of differential di-jet cross-sections, 
and using both the inclusive k±_ and a cone algorithm for the study of jet structure. 

Jet shapes, ^(r), have been studied in two separate samples: xz > 0.75, which is dom- 
inated by direct processes and hence by quark-initiated jets, and x^ < 0.75, dominated by 
resolved events and therefore by gluon-initiated jets. As expected from QCD the jets in 
the first sample are significantly more collimated than for x^ < 0.75. Jets in both sam- 
ples become more collimated with increasing transverse energy, but show no significant 
dependence on the jet pseudo-rapidity. Jets defined by the cone algorithm are substan- 
tially broader than those defined by the k± algorithm at low E^ . However the difference 
decreases with increasing E^. The shape of fc^-jets is well described by PYTHIA and 
PHOJET. The jet shapes obtained for the cone-jets are somewhat broader in PYTHIA 
and PHOJET than in the data. 

Inclusive differential di-jet cross-sections have been measured as a function of |cosO*|, 
E^, \rj^ et \ and lA?^ 6 *! and, for the first time, as a function of x 7 in several bins of E^ . 
Different regions of the £+-2~-space are explored to separate experimentally direct from 
resolved interactions and to study and control the influence of an underlying event. By 
measuring the cross-sections for events in which either x+ or x~ is smaller than 0.75 we have 
isolated a region of phase space in which resolved photon processes dominate, and which 
at the same time is much less sensitive to multiple parton interactions. By performing the 
measurement also for x^ < 0.75, observables are made available which are sensitive to the 
amount of multiple parton interactions added in the prediction, and which can be used to 
study these effects in detail. 

A strong rise with increasing |cosO*| is observed for the differential di-jet cross-section 
for x^ < 0.75, as expected from QCD for a sample with a significant contribution from 
spin-1 gluon exchange. The flatter distribution for direct events is also in good agreement 
with the QCD calculation. 
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The differential di-jet cross-sections as a function of E 3 ^ , \if et \ : |A?/j ct | and x 1 are in 
good agreement with the next-to-leading order perturbative QCD calculation except for 
x^ < 0.75, where the calculation is too low. As this calculation does not include a model 
for the underlying event that is expected to be largest in this region, it is interesting to note 
that the discrepancy is of similar size to the contribution of multiple parton interactions 
to the PYTHIA prediction. The sensitivity of the results presented to the gluon density in 
the photon is clearly visible in NLO QCD predictions using different parton distribution 
functions, but is compensated to some extent by anti-correlated differences in the respective 
quark-densities. A global analysis using additional data sets to simultaneously constrain 
the quark densities hence promises to yield the highest sensitivity to the gluon density in 
the photon. 

The measurements carried out for events in which only either x^ or x~ is smaller than 
0.75 are a unique data set. While this region is almost insensitive to multiple parton 
interactions, the fraction of events at small x 7 is still sizable, which indicates a significant 
contribution of resolved processes and hence a good sensitivity to the hadronic structure of 
the photon. The good agreement of data and theory in this region in particular confirms 
that perturbative QCD in next-to-leading order is able to describe correctly the inclusive 
production of di-jets in photon-photon collisions. 
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cos0* 



d|cos6* 



[Pb] 



' stat 



[Pb] 



'sys 



[Pb] 



x± > 0.75 


0.0 


0.1 


5.04 




0.49 


0.32 


0.1 


0.2 


4.96 




0.49 


0.34 


0.2 


0.3 


5.22 




0.49 


0.32 


0.3 


0.4 


5.35 




0.48 


0.38 


0.4 


0.5 


4.90 




0.45 


0.32 


0.5 


0.6 


6.73 




0.56 


0.47 


0.6 


0.7 


7.09 




0.54 


0.43 


0.7 


0.8 


8.42 




0.61 


0.51 


x± < 0.75 


0.0 


0.1 


1.50 




0.31 


0.24 


0.1 


0.2 


2.63 




0.46 


0.36 


0.2 


0.3 


2.91 




0.50 


0.41 


0.3 


0.4 


3.73 




0.62 


0.43 


0.4 


0.5 


4.36 




0.61 


0.62 


0.5 


0.6 


6.52 




0.80 


0.74 


0.6 


0.7 


9.90 




0.92 


1.07 


0.7 


0.8 


13.05 




0.96 


1.02 



Table 1: The di-jet cross-section as a function of |cos©*| for the two regions in x+-x~- 
space indicated in the table. The total uncertainty for each bin is the quadratic sum of the 
statistical and systematic uncertainty given in the table. 
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[GeV] 


ST [Pb/Gev] 


5 s tat [pb/GeV] 


5 sys [pb/GeV] 


full - x 7 - range 


5.00 6.54 

6.54 8.55 

8.55 11.18 
11.18 14.62 
14.62 19.12 
19.12 25.00 


11.54 
4.16 
1.45 
0.543 
0.176 
0.0607 


0.36 
0.16 
0.07 
0.036 
0.020 
0.0117 


0.59 
0.25 
0.09 
0.044 
0.020 
0.0087 


x+ or x^ < 0.75 


5.00 6.54 

6.54 8.55 

8.55 11.18 
11.18 14.62 
14.62 19.12 


3.56 
1.32 
0.475 
0.190 
0.0518 


0.21 
0.09 
0.045 
0.025 
0.0127 


0.22 
0.10 
0.037 
0.024 
0.0073 


x$ < 0.75 


5.00 6.54 

6.54 8.55 

8.55 11.18 
11.18 14.62 


5.94 
1.84 
0.461 
0.0781 


0.41 
0.16 
0.060 
0.0192 


0.66 
0.23 
0.071 
0.0130 



Table 2: The di-jet cross-section as a function of the mean transverse energy of the 
di-jet system, for the three regions in x+-x~ -space indicated in the table. 
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[pb] 


8 stat [pb] 


5 sys [pb] 




full x+ - x^ - 


range 




5 GeV < < 7 GeV 


0.000 


0.125 


47.29 


2.51 


4.41 


0.125 


0.250 


46.29 


2.73 


4.69 


0.250 


0.375 


33.66 


2.16 


3.37 


0.375 


0.500 


25.47 


1.93 


2.25 


0.500 


0.625 


25.70 


1.86 


2.16 


0.625 


0.750 


38.91 


2.47 


2.70 


0.750 


0.875 


79.59 


4.07 


8.37 


0.875 


1.000 


51.97 


5.12 


5.49 


7 GeV < < 11 GeV 


0.000 


0.125 


12.41 


1.37 


1.70 


0.125 


0.250 


16.78 


1.55 


2.05 


0.250 


0.375 


12.64 


1.34 


1.47 


0.375 


0.500 


10.36 


1.16 


0.86 


0.500 


0.625 


11.67 


1.24 


0.95 


0.625 


0.750 


14.80 


1.53 


1.05 


0.750 


0.875 


38.84 


2.72 


3.99 


0.875 


1.000 


43.07 


3.00 


4.75 


11 GeV < E^ < 25 GeV 


0.000 


0.125 


1.51 


0.53 


0.39 


0.125 


0.250 


2.83 


0.52 


0.64 


0.250 


0.375 


2.78 


0.54 


0.43 


0.375 


0.500 


2.30 


0.48 


0.36 


0.500 


0.625 


2.97 


0.57 


0.41 


0.625 


0.750 


4.43 


0.66 


0.55 


0.750 


0.875 


11.04 


1.23 


0.94 


0.875 


1.000 


26.81 


1.85 


2.09 



Table 3: The di-jet cross-section as a function of x 1 for the regions of the mean transverse 
energy E^ of the di-jet system and the region in x+-x~ -space indicated in the table. 
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[pb] 


5 stat [pb] 


5 sys [pb] 


x+ or x^ < 0.75 


5 GeV < < 7 GeV 


0.000 


0.125 


12.47 


0.96 


1.18 


0.125 


0.250 


8.87 


0.74 


1.16 


0.250 


0.375 


7.00 


0.69 


0.83 


0.375 


0.500 


5.12 


0.58 


0.65 


0.500 


0.625 


5.99 


0.69 


0.56 


0.625 


0.750 


12.36 


1.27 


0.98 


0.750 


0.875 


38.50 


1.92 


4.19 


0.875 


1.000 


11.80 


1.72 


1.30 


7 GeV < < 11 GeV 


0.000 


0.125 


5.44 


0.73 


0.62 


0.125 


0.250 


4.84 


0.52 


0.48 


0.250 


0.375 


3.45 


0.44 


0.31 


0.375 


0.500 


2.91 


0.41 


0.21 


0.500 


0.625 


3.28 


0.48 


0.25 


0.625 


0.750 


5.38 


0.76 


0.39 


0.750 


0.875 


15.03 


1.17 


1.94 


0.875 


1.000 


8.63 


1.19 


1.31 



Table 4-' The di-jet cross-section as a function of x~ t for the regions of the mean transverse 
energy E^ 1 of the di-jet system and the region in x^-x^ -space indicated in the table. 
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tpb] 


5 stat [pb] 


5 sys [pb] 


x± < 0.75 


5 GeV < < 7 GeV 


0.000 


0.125 


42.37 


2.94 


4.71 


0.125 


0.250 


37.66 


2.44 


4.32 


0.250 


0.375 


26.87 


1.99 


2.66 


0.375 


0.500 


19.63 


1.88 


1.49 


0.500 


0.625 


17.84 


1.55 


1.10 


0.625 


0.750 


26.71 


3.86 


1.73 


7 GeV < < 11 GeV 


0.000 


0.125 


8.81 


1.08 


1.54 


0.125 


0.250 


11.00 


1.00 


1.72 


0.250 


0.375 


10.30 


0.96 


1.24 


0.375 


0.500 


8.10 


0.81 


0.76 


0.500 


0.625 


7.49 


0.94 


0.64 


0.625 


0.750 


8.51 


1.58 


0.72 



Table 5: The di-jet cross-section as a function of x 7 for the regions of the mean transverse 
energy E^ of the di-jet system and the region in x+-x~ -space indicated in the table. 
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logi (x 




dCTdijct [pbl 


8 stat [pb] 


5 sys [pb] 


7) 


dlogio(a; 7 ) L-f J 




5 GeV < 


< 7 GeV 








full x+ - x 1 


- ranj 


?e 




-1.65 


-1.32 


4.00 






0.42 


0.42 


-1.32 


-0.99 


9.54 






0.76 


1.07 


-0.99 


-0.66 


16.55 






1.12 


1.83 


-0.66 


-0.33 


27.51 






1.39 


2.35 


-0.33 


0.00 


76.87 






1.69 


5.05 






x+ or x~ 


< 0.75 






-1.65 


-1.32 


1.29 






0.18 


0.18 


-1.32 


-0.99 


2.38 






0.24 


0.32 


-0.99 


-0.66 


3.64 






0.30 


0.44 


-0.66 


-0.33 


4.87 






0.38 


0.47 


-0.33 


0.00 


25.87 






0.99 


2.16 


x$ < 0.75 


-1.65 


-1.32 


3.08 






0.44 


0.41 


-1.32 


-0.99 


7.81 






0.63 


1.02 


-0.99 


-0.66 


13.95 






0.92 


1.55 


-0.66 


-0.33 


19.39 






1.13 


1.78 


-0.33 


0.00 


18.01 






1.53 


1.36 



Table 6: The di-jet cross-section as a function of logio(x 7 ) for the region of the mean 
transverse energy of the di-jet system and the regions in xt f -x~ -space indicated in the 
table. 
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d^dijct r ii 


5 stat [pb] 


5 sys [pb] 






x+ or x~ 


< 0.75 




0.0 


0.5 


6.68 


0.27 


0.38 


0.5 


1.0 


5.51 


0.24 


0.31 


1.0 


1.5 


3.98 


0.20 


0.23 


1.5 


2.0 


2.34 


0.15 


0.15 


2.0 


2.5 


1.22 


0.11 


0.11 


2.5 


3.0 


0.631 


0.081 


0.074 


3.0 


3.5 


0.251 


0.056 


0.035 


x± < 0.75 


0.0 


0.5 


7.47 


0.36 


0.66 


0.5 


1.0 


6.89 


0.35 


0.60 


1.0 


1.5 


5.42 


0.29 


0.46 


1.5 


2.0 


3.93 


0.24 


0.30 


2.0 


2.5 


2.44 


0.19 


0.30 


2.5 


3.0 


1.43 


0.16 


0.24 


3.0 


3.5 


0.618 


0.124 


0.133 



Table 7: The di-jet cross-section as a function of |A?7j ct | for the two leading jets in , 
for the two regions in x~t-x~Z -space indicated in the table. 
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^cntr 






S stat [Pb] 


5 sys [pb] 


x+ or x 7 < 0.75 


0.00 0.25 
0.25 0.50 
0.50 0.75 
0.75 1.00 
1.00 1.25 
1.25 1.50 
1.50 1.75 


11.48 

9.25 

7.19 

5.30 

3.49 

2.60 

1.49 


0.50 
0.44 
0.38 
0.32 
0.26 
0.24 
0.21 


0.65 
0.55 
0.46 
0.31 
0.23 
0.25 
0.17 


x± < 0.75 


0.00 0.25 
0.25 0.50 
0.50 0.75 
0.75 1.00 
1.00 1.25 
1.25 1.50 
1.50 1.75 


14.93 

13.38 

10.83 

7.17 

5.59 

3.20 

2.01 


0.67 
0.65 
0.60 
0.47 
0.46 
0.41 
0.41 


1.24 
1.17 
0.96 
0.60 
0.54 
0.41 
0.32 



Table 8: The di-jet cross-section as a function of l^ntrl f or ^ e ^ wo regions in x^-x^ -space 
indicated in the table. 
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1 j ct 1 
l 7 7fwdl 




tfstat [pb] 


5 sys [pb] 


x+ or x 7 < 0.75 


0.00 




0.25 


0.906 


0.137 


0.065 


0.25 




0.50 


2.68 


0.23 


0.18 


0.50 




0.75 


4.30 


0.29 


0.26 


0.75 




1.00 


5.65 


0.33 


0.33 


1.00 




1.25 


6.18 


0.34 


0.33 


1.25 




1.50 


6.86 


0.36 


0.34 


1.50 




1.75 


7.27 


0.40 


0.48 


1.75 




2.00 


7.03 


0.44 


0.55 


x± < 0.75 


0.00 




0.25 


1.03 


0.17 


0.16 


0.25 




0.50 


4.11 


0.38 


0.54 


0.50 




0.75 


5.64 


0.40 


0.60 


0.75 




1.00 


6.71 


0.41 


0.51 


1.00 




1.25 


8.40 


0.47 


0.74 


1.25 




1.50 


10.72 


0.62 


0.94 


1.50 




1.75 


11.30 


0.68 


1.18 


1.75 




2.00 


9.86 


0.72 


0.96 



Table 9: The di-jet cross-section as a function of \rf^ d \ for the two regions in x+-x^ -space 
indicated in the table. 
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Figure 1: (a)-(c): The relative contribution of direct, single-resolved, and double-resolved 
processes according to PYTHIA at the hadron level for the cross-sections as a function 
of Xry for the full x+ -x~ -range (see Section^. In (d) the regions in x+-x~ -space used in 
addition to the full x+ -x~ -range (referred to as (A)) are illustrated: (B) both x+ and x~ 
larger than 0.75 (x^ > 0.75), (C) either x+ or x~ smaller than 0.75 (x+ or x~ < 0.75), 
(D) both and x~ smaller than 0.75 (x^ < 0.75). 



30 




Figure 2: Comparison of event quantities for uncorrected data with the simulation for 
the di-jet sample including contributions from background processes, (a) shows the sum of 
energy measured in the electromagnetic and hadronic calorimeter, (b) the number of tracks 
in the event, (c) is the sum of energy in the FD detector scaled by the beam energy, and 
(d) the invariant mass of the jet with the highest in the event and the four-vector 
calculated from all objects in the opposite hemisphere as seen from this jet. The statistical 
error is shown where larger than the marker size. The label 7*7 stands for simulated 
photon-photon collision events in which one of the photons has a virtuality larger than 
4-5 GeV 2 as discussed in Section^ 
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Figure 3: Jet profiles: The Et-JIow normalised to the number of jets Nj et as a function 
of the distance from the jet axis in <fi o,nd fj. Jets are selected from the range 10 < Ei£ < 
25 GeV. The statistical error is shown where larger than the marker size. The data are 
compared to a mixture of signal (PHOJET or PYTHIA) and background MC simulation. 
The background MC simulations used are the same as in Figure® 
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Figure 4- The energy-flow as a function off} where the area around the two leading jets 
in the event has been excluded in a cone of radius 1.3 in rj-(j)-space. The distributions are 
normalised to the number of di-jet events N. The statistical error is shown where larger 
than the marker size. The regions in xt f -x~ -space considered and the MC simulations used 
are the same as in Figured 
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Figure 5: The uncorrected x 7 distributions in the data compared to the sum of signal and 
background processes in the simulation. The statistical error is shown where larger than 
the marker size. The MC simulations used are the same as in Figure® 
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Figure 6: Hadronisation corrections estimated by PYTHIA and HERWIG for (a) Ei£ and 
(b)-(d) Xj for the regions of E^ given in the figure. In all cases the full x+-x~ -space is 
considered. 
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Figure 7: The jet shape, ^(r), for the two regions of x^ -x~ -space indicated in the figure (a), 
and \I/(r) for < 0.75 compared to the predictions of the LO MC generators PHOJET 
and PYTHIA (b). Figures (c) and (d) show the value of ^(r = 0.4) as a function of 
the transverse energy and pseudo-rapidity of the jet respectively, compared to the PYTHIA 
prediction. The total of statistical and systematic uncertainties added in quadrature is 
shown where larger than the marker size. The inner error bars show the statistical errors. 
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Figure 8: The value of the jet shape \P(r) at r = 0.4 as a function of the jet transverse 
energy for < 0.75 (a) and > 0.75 (b), and as a function of the jet pseudo-rapidity 
forx^ < 0.75 (c) and x^ > 0.75 (d). In each figure the results obtained using the inclusive 
k± and the cone jet algorithm are shown and compared to the PYTHIA prediction. The 
total of statistical and systematic uncertainties added in quadrature is shown where larger 
than the marker size. The inner error bars show the statistical errors. 
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Figure 9: The di-jet cross-section as a function of |cos6*| for the two regions in x+-x~- 
space indicated in the figure. The total of statistical and systematic uncertainties added 
in quadrature is shown where larger than the marker size. The inner error bars show 
the statistical errors. The numerical values are given in TableU\ A perturbative NLO 
QCD prediction ' t 29j using the GRVHO parton densities is compared to the data after 
hadronisation corrections have been applied to the calculation. The shaded band indicates 
the theoretical uncertainty estimated by the quadratic sum of two contributions: variation of 
the renormalisation scale by factors of 0. 5 and 2 and the difference between using HERWIG 
or PYTHIA in estimating the hadronisation corrections. 
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Figure 10: The di-jet cross-section as a function of the mean transverse energy of the 
di-jet system, for the three regions in x+-x~ -space given in the figure. The factor f is used 
to separate the three measurements in the figure more clearly. The total of statistical and 
systematic uncertainties added in quadrature is shown where larger than the marker size. 
The inner error bars show the statistical errors. The numerical values are given in Tabled 
The prediction of the LO program PYTHIA using the parton distribution function SaS ID 
is compared to the data. The NLO calculation is the same as in Figure^ 
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Figure 11: The di-jet cross-section as a function of x 1 and logio (x 7 ) for the regions of 
the mean transverse energy of the di-jet system indicated in the figures. The total 
of statistical and systematic uncertainties added in quadrature is shown where larger than 
the marker size. The inner error bars show the statistical errors. The numerical values 
are given in Tables® and® The prediction of the LO MC generator PYTHIA using the 
parton distribution function SaS ID is compared to the data. The shaded histogram at the 
bottom of each plot indicates the MIA contribution to the PYTHIA prediction. The NLO 
calculation is the same as in Figure® 
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Figure 12: The prediction of NLO QCD using different parton distributions for the photons. 
The di-jet cross-section as a function ofx 1 for 5 GeV < < 7 GeV (upper left plot in 
Figure [77)) is shown. In (a) the full cross-section is shown after hadronisation corrections 
have been applied, while (b), (c) and (d) show the gg, gq and qq contributions to this 
cross-section without hadronisation corrections. The NLO calculation is the same as in 
Figure fji 
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Figure 13: The di-jet cross-section as a function of x^ and logio (x 7 ) for the regions of 
the mean transverse energy of the di-jet system indicated in the figures and for x+ or 
x~ < 0.75 . The total of statistical and systematic uncertainties added in quadrature is 
shown where larger than the marker size. The inner error bars show the statistical errors. 
The numerical values are given in Tables^and\$ The NLO calculation and MC simulation 
are the same as in Figure 1771 
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Figure 14.: The di-jet cross-section as a function of x 7 and logio (%) / or regions of the 
mean transverse energy of the di-jet system indicated in the figures. For these cross- 
sections < 0.75 is required. The total of statistical and systematic uncertainties added 
in quadrature is shown where larger than the marker size. The inner error bars show the 
statistical errors. The numerical values are given in Tables^ and\$ The NLO calculation 
and MC simulation are the same as in Figure [771 
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Figure 15: The di-jet cross-section as a function of \Arf et \ for the two leading jets in 
and separately for the central and the forward jet, for x+ or x~ < 0.75. The total of 
statistical and systematic uncertainties added in quadrature is shown where larger than the 
marker size. The inner error bars show the statistical errors. The numerical values are 
given in Tables\7\\$ and\Q The NLO calculation and MC simulation are the same as in 
Figure [771 
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Figure 16: The di-jet cross-section as a function of |A^ et | for the two leading jets in E^ t 
and separately for the central and the forward jet. For these cross-sections < 0.75 is 
required. The total of statistical and systematic uncertainties added in quadrature is shown 
where larger than the marker size. The inner error bars show the statistical errors. The 
numerical values are given in Tables^7\\$ and\Q The NLO calculation and MC simulation 
are the same as in Figure [771 
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